Atmospheric neutrino observations and flavor changing interactions 
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Flavor changing (FC) neutrino-matter interactions can account for the zenith-angle dependent 
deficit of atmospheric neutrinos observed in the SuperKamiokande experiment, without directly in- 
voking neither neutrino mass, nor mixing. We find that FC //^-matter interactions provide a good 
fit to the observed zenith angle distributions, comparable in quality to the neutrino oscillation hy- 
pothesis. Fhe required FC interactions arise naturally in many attractive extensions of the Standard 
Model. 



Neutrinos produced as decay products in hadronic 
showers from cosmic ray collisions with nuclei in the up- 
per atmosphere [jlj have been observed by several detec- 
tors . Although the absolute fluxes of atmospheric 
neutrinos are largely uncertain, the expected ratio (/i/e) 
of the muon neutrino flux (y^ + D^) over the electron 
neutrino flux (y e + v e ) is robust, since it largely cancels 
out the uncertainties associated with the absolute flux. 
In fact, this ratio has been calculated JjJ with an uncer- 
tainty of less than 5% over energies varying from 0.1 GeV 
to 100 GeV. In this resides our confidence on the long- 
standing atmospheric neutrino anomaly. 

Although the first iron-calorimeter detectors in 
Frejus @] and NUSEX g reported a value of the dou- 
ble ratio, K(fx/e) = (/V e )data/(/V e )MCj consistent with 
one, all the water Cherenkov detectors Kamiokande [Q, 
1MB ]{| and SuperKamiokande |J have measured R(/x/e) 
significantly smaller than one. Moreover, not long 
ago, the Soudan-2 Collaboration, also using an iron- 
calorimeter, reported a small value of R(/i/e) j7j], showing 
that the so-called atmospheric neutrino anomaly was not 
a feature of water Cherenkov detectors. 

Recent SuperKamiokande high statistics observa- 
tions H indicate that the deficit in the total ratio R(/x/e) 
is due to the number of neutrinos arriving in the detec- 
tor at large zenith angles. Although e-like events do not 
present any compelling evidence of a zenith-angle depen- 



dence, the /i-like event rates are substantially suppressed 
at large zenith angles. 

The — > v r [||H as well as the — > v s |||| oscilla- 
tion hypothesis provide an appealing explanation for this 
smaller-than-expected ratio, as they are simple and well- 
motivated theoretically. This led the SuperKamiokande 
Collaboration to conclude that their data provide good 
evidence for neutrino oscillations and neutrino masses. 

In this letter we give an alternative explanation of 
the atmospheric neutrino data in terms of FC neutrino- 
matter interactions @ 0. We show that even if neu- 
trinos have vanishing masses and/or the vacuum mixing 
angle is negligible, FC neutrino matter interactions can 
still explain the SuperKamiokande data. 

There are attractive theories beyond the SM where 
neutrinos are naturally massless Jl5[ as a result of a pro- 
tecting symmetry, such as B-L in the case of supersym- 
metric SU(5) models Q and the model proposed in |l7]] 
or chiral symmetry in theories with extended gauge struc- 
ture such as SU(3) c (g>SU{3) L (g>U(l) N (331) models [g. 
The simplest example of this mechanism was first noted 
in an SU(2) x U(l) model proposed in Ref. |l^,[l7| where 
a singlet Dirac lepton at the TeV scale is added sequen- 
tially to the SM in such a way that neutrinos remain 
massless due to an imposed B-L symmetry. The flavor 
mixing amongst the massless neutrinos in the leptonic 
charged current can not be rotated away despite neutri- 
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nos being massless. 

In general, models of neutrino mass are a natural 
source of FC neutrino-matter interactions. Seesaw-type 
models of neutrino mass have non-diagonal neutral cur- 
rent couplings of the Z to mass eigenstate neutrinos ]l9| ] 
that may lead to new FC neutrino-matter interactions. 
Models of radiative generation of neutrino mass [^0| typ- 
ically contain additional FC neutrino-matter interactions 
from scalar exchanges. An example of this class are the 
supersymmetric models with broken i?-parity pl|] . Mod- 
els with extended gauge structure, such as E(6) models 
jn|, may also lead to FC neutrino-matter interactions. 

Here we focus on a massless neutrino conversion sce- 
nario as an explanation of the atmospheric neutrino data 
based on FC neutrino-matter interactions which induce 
— > v T transition. Our results can be extended to 
the massive neutrino case. However, we stress that the 
present atmospheric neutrino data does not necessarily 
provide evidence for neutrino mass. Moreover, the exis- 
tence of attractive theories where FC neutrino-matter in- 
teractions do not imply neutrino mass, makes this possi- 
bility especially elegant. From a phcnomcnological point 
of view FC interactions of neutrinos can induce flavor 
transitions when neutrinos travel through matte r |10[| ir- 
respective of neutrino mass. In both massive |22f and 
massless Jl^,|l3|] cases conversions can be resonant, how- 
ever the properties of the conversion are totally differ- 
ent. Massless neutrino conversions would be energy- 
independent and would affect atmospheric neutrinos as 
well as anti-neutrinos, converting — > v T together with 

—* D T . Remarkably the present double ratio data does 
not show a significant energy dependence |6|,[23fl- 

The presence of FC neutrino-matter interactions im- 
plies a non-trivial structure of the neutrino evolution 
Hamiltonian in matter. The evolution equations describ- 
ing the Ufj, — ► v T transitions in matter are given as fll2| , |i"3| : 
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where, v a = v a (r), a=/x, t are the probability amplitudes 
to find these neutrinos at a distance r from their creation 
position, V2 Gpnf{r)e v is the + / — > v T + / forward 
scattering amplitude and y/2 GFnf(r)e' v is the difference 
between the v T — f and — / elastic forward scattering 
amplitudes, with n/(r) being the number density of the 
fermions which induce such processes. 

The use of the FC i/ e -matter interactions was previ- 
ously suggested in connection with the solar neutrino 
problem Recently an attempt was made to 

extend this idea in order to account also for the atmo- 
spheric neutrino data, but the fit obtained in this paper 
is not as good as our atmospheric neutrino fit or the solar 
fit in Ref. KM. Moreover, some of the results in Table 
III of Ref. |24| seem inconsistent. On the other hand 
Ref. E3] includes exotic flavor-conserving v T interactions 



plus neutrino masses in order to account for the atmo- 
spheric and LSND data, but without a detailed fit. We 
have decided to postpone the detailed analysis (within 
the present scenario) of the LSND, as well as of the solar 
neutrino data, for a future work. Here we show that FC 
i^-matter interactions can explain the atmospheric neu- 
trino zenith angle anomaly, without introducing neutrino 
masses and/or mixing. 

For our phenomenological approach let us simply as- 
sume the existence of a tree-level process v a + f + / 
with amplitude proportional to g a fg0f/4m 2 , where a 
and (3 are flavor indices, / stands for the interacting ele- 
mentary fermion (charged lepton, d-like or w-like quark) 
and g a f is the coupling involved in the vertex where a v a 
interacts with / through a scalar or vector boson of mass 
to. The evolution equations which describe the — > v T 
transitions in matter may be written generically as in 
Eq. (nl), where: 



4to 2 V2Gf 
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Note also that, in the absence of neutrino mass as we are 
assuming, anti-neutrino transitions — > v r are governed 
by precisely the same evolution matrix in Eq. ([l]). We 
have calculated the transition probabilities of — > v T 
(z? M — * D T ) P(e v ,e' u ), (P(e p ,£p)) as a function of the 
zenith angle by numerically solving the evolution equa- 
tion using the density distribution in and a realistic 
chemical composition with proton/neutron ratio 0.497 
in the mantle and 0.468 in the core p^| . For the sake 
of simplicity we have assumed that e = e p = e„ and 
e' = e' B = ej,, so there are only two free parameters in 
the analysis. We have used these probabilities to com- 
pute the theoretically expected number of fi- and e-like 
events (N^ and N e ) as a function of the two parameters, 
e and e', for each of the 5 zenith angle bins both for the 
sub-GeV and for the multi-GeV SuperKamiokande data. 
Following Rcfs. p^8|, we fit separately the \x- and e-like 
events taking into account the correlation of errors. The 
calculated numbers of events, and 7V e , as functions 
of zenith angle, have been compared with the 535-days 
SuperKamiokande data sample in order to determine the 
allowed regions of e and e' from a \ 2 fit. We set our 
normalization assuming that the relevant neutrino inter- 
action in the Earth is only with down- type quarks. Any 
other scenario can be obtained from our results by rescal- 
ing the e-parameters. 

In Fig.([l]) we present a contour plot of the regions al- 
lowed by the SuperKamiokande data. The contour plots 
(a), (b) and (c) correspond to the regions allowed by the 
sub-GeV, multi-GeV and combined SuperKamiokande 
data, respectively. These contours are determined by the 
conditions \ 2 = Xmin + ^X 2 where A\ 2 = 4.6, 6.0, 9.2 
for 90, 95 and 99 % C. L., respectively. 



In the parameter region we have considered, i.e., e and 
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FIG. 1. Allowed region for e and e for SuperKamiokande 
(a) sub-GeV (b) multi-GeV and (c) combined events in the 
massless-neutrino scenario. The best fit points for each case 
is indicated by the crosses. 



e' in the interval [0.01, 1.0], we found that Xmin = 6-3 
and 6.4 for the sub-GeV and multi-GeV samples (8 d. o. f. 
corresponding to 10 data points minus two free parame- 
ters). These minima are obtained for (e, e') = (0.08, 0.01) 
and (0.68, 0.36), respectively, as indicated by the crosses 
in Fig.(|l|). For the combined case, Xmin — 14.7 (18 d.o.f) 
for (e, e') = (0.99, 0.02). In the combined case the local 
best fit point { X 2 = 16.9 for (e, e') = (0.08, 0.07) ) in 
the "island" determined by the 90 % C. L. curve is also 
indicated by a filled circle. This point is interesting be- 
cause it still gives a good fit to the data with a relatively 
small value for the FC parameter e. We find also that 
the x 2 is relatively flat along the e' axis around the best 
fit point. The allowed regions can be qualitatively under- 
stood in the approximation of constant matter density. 
The conversion probability in this case is: 



ie 2 



4e 2 



r'2 



sin 2 (-r]L), 



t' 2 ^/2GFnf . For nj — rid ?s 3n e 
< e, the oscillation length in matter is given by: 
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and 
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From Eq. (|) one can see that in order to have a large 
transition probability one must be in the region e' < e 
and rj < tt/Rq. This last condition leads to a lower 
bound on e. The island in Fig. l.(b) corresponds to 

In Fig.(||) we give the expected zenith angle distribu- 
tion of /x-like sub-GeV events (a) and multi-GeV events 
(b) evaluated with our Monte Carlo program for the best 
fit points determined above. Our results clearly indicate 
an excellent fit for the /z-like events showing that they 
are highly depleted at cos0 = — 1 with respect to the SM 



prediction. Note that, except for the assumption that the 
FC ^-matter interaction involves d-quarks, our result is 
quite general, since we have not explicitly considered any 
particular model as the origin of the FC neutrino-matter 
interaction. Note that e-like events are not affected by 
— > v T transition. 

What can we say about the required strength of the 
neutrino-matter interaction in order to obtain a good 
fit of the observed data? From our results and Eq. (|^) 
we see that for masses m « 200 GeV we need at least 
g T f ■ g^f ^0.1 for the the mixing term e. Similarly our 
best fit e' value implies \g T f\ 2 — 1 5^/1 2 ~ 0.1. While these 
values are relatively large, they are both weak-strength 
couplings. Moreover they are consistent with present ex- 
perimental bounds, for example from universality of the 
weak interaction which is manifestly violated by Eq. (Q). 

For the purpose of illustrating this explicitly let us con- 
sider for the moment the supersymmetric model with 
broken i?-parity as a way to parameterize the FC 
neutrino-matter interaction. In this case the FC v^- 
matter interactions are mediated by a scalar down-type 
quarks, dj, so that we need only to check the couplings 
where a d-quark and a fi- or r-neutrino is involved, i.e 
9id ~ Kjii « = 2,3. \lj k are the coupling constants in 
the broken i?-parity superpotential X' i , k LiQjD^ where 



L, Q and D are standard superfields, and A^/2Gf^ 

3 




-0.5 0.0 0.5 
COS© 



0.5 0.0 0.5 
COS© 



FIG. 2. Best-fit zenith angle distributions in the mass- 
less-neutrino FC scenario (thin-solid, dotted lines) ver- 
sus no-oscillation hypothesis (thick-solid line). The Su- 
perKamiokande data are indicated by the crosses. 

Constraints on the magnitude of such FC interactions 
in broken i?-parity models have been given in p^j. The 
precision tests imply that individually these couplings are 
not severely constrained. The most stringent limit to the 
values of the relevant FC quantities comes from limits on 



the FC tau decay BR(t 



(j,-) < 6.3 x 10~ 6 EC 
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which implies that | V]A^ 1 A^ 1 (100 GeV/m^) 2 ! < 3.1 x 

3 



10~ 3 . Although a certain degree of fine-tuning is needed 
in order to verify this constraint we find that the required 
strength of FC ^-matter interaction is consistent with 
all present data. One should note that broken i?-parity 
supersymmctric models typically lead to neutrino masses 
which could be large. However one may suppress them 
via fine-tuning or some additional symmetry. 

In summary, we have demonstrated that flavor chang- 
ing i^-matter interactions can account for the zenith- 
dependent deficit of atmospheric neutrinos observed in 
the SuperKamiokande experiment, without directly in- 
voking neutrino masses and mixing. It provides a fit of 
the observations which is significantly better than the no- 
oscillation hypothesis and of similar quality as the usual 
— ► v r oscillation hypothesis. The required FC inter- 
action can arise in many attractive extensions of the SM 
and is consistent with all present constraints. 

The above FC mechanism can also be tested at future 
Long Baseline experiments. From Eq. (^), using n e ~ 2 
mol/cc, e ~ 1 (0.1) and e' < e, for the planned K2K 
experiment |51j one gets Piv^ — > v T ) ~ 0.4 (0.004) while 
for MINOS |f one finds P(i^ -> v T ) ~ 0.9 (0.2). 

The existence of a massless neutrino explanation of the 
atmospheric neutrino anomaly may play an important 
theoretical role in model-building, especially if one wants 
to account for all other hints for non-standard neutrino 
properties, namely the solar neutrino data, the LSND 
result, and the possible role of neutrinos as dark matter 
in the Universe |T^ |. 
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